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Classification of wood surface texture based on Gauss-MRF Model 


WANG Ke-qi, BAI Xue-bing 
Northeast Forestry University, Harbin 150040, P. R. China 


Abstract: The basal theory of Gauss-MRF is expounded and 2-5 order Gauss-MRF models are established. Parameters of the 2-5 order 
Gauss-MRF models for 300 wood samples’ surface texture are also estimated by using LMS. The data analysis shows that: 1) different tex¬ 
ture parameters have a clear scattered distribution, 2) the main direction of texture is the direction represented by the maximum parameter of 
Gauss-MRF parameters, and 3) for those samples having the same main direction, the finer the texture is, the greater the corresponding pa¬ 
rameter is, and the smaller the other parameters are; and the higher the order of Gauss-MRF is, the more clearly the texture is described. On 
the condition of the second order Gauss-MRF model, parameter B1, B2 of tangential texture are smaller than that of radial texture, while B3 
and B4 of tangential texture are greater than that of radial texture. According to the value of separated criterion, the parameter of the fifth 
order Gauss-MRF is used as feature vector for Hamming neural network classification. As a result, the ratio of correctness reaches 88%. 
Keywords: Wood surface texture; Gauss-MRF; Feature parameter; Parameter estimation, Separation judgment; Classification 
C1C number: S781.1 Document code: A Article ID: 1007-662X(2006)01-0057-05 


Introduction 

The natural texture has a complicated structure. Presently 
there has not been a standard to describe its feature. How to de¬ 
scribe and classify natural textures is an important issue in the 
field of computer vision (Li and Liu 1998). The description of 
texture feature is divided into statistical method, structural 
method, and model method (Sheng 2000) It has been proved by 
practice that the statistical method represented by Gray Level 
Co-occurrence Matrix is suitable for dealing with the natural 
texture, but its feature parameter is unable to clearly indicate the 
geometrical feature of the texture and the computation work is 
too large. The random model method represented by Markov 
random field has a unique advantage in this aspect because the 
parameters produced by the method can describe the gathering 
characters of the texture in different directions and different pat¬ 
terns (Gong 1999), which is in good accordance with the sensory 
organ of human beings. 

Classification of texture remains with pattern-recognition. 
Among those texture feature parameters, there exists both a de¬ 
pendent relation and a non-linear relation. Meanwhile, it involves 
an estimation issue of priori probability density. Therefore it is 
difficult to classify the texture with traditional Bayes decision 
theory (Bian et al. 2000). Being similar to the human brain in 
processing complicated non-linear information, the neural net¬ 
work brings about a vital force for pattern-recognition field, of 
which the most noticeable ones are BP neural network and com¬ 
petition neural network. 

A method with Gauss-Markov Random Field to set up feature 
parameters of texture is put forward in this paper, and the natural 
texture of wood is classified with the competition neural net¬ 
work. 
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Gauss-Markov Random Field 

Markov Random Field (abbreviated as MRF) is the mathe¬ 
matics theory of random process description put forward by D. 
Geman and S. Geman in 1984. Its definition is as follows: 

Let S = (Si, S 2 ,...,S^) be the position set in a 2D space, Abe 
the random variable. The random event (xi, x 2 , ... , x N ) can be 
expressed as (X=a>) , co£f2, Q is the value space of random 
variable; C, is the neighborhood system of the 2D space. All only 
and, CO e Q , p(X = CO) > 0 , we have 

p(x s \x r ,r* s) = p(x s \x r ,re |) (1) 

For all S, (O holds, then {X} is a Markov random field. A 
2D picture can be regarded as a Markov random field, which 
means that the gray level probability of any picture element in a 
picture field is decided by the gray level of picture element in its 
neighborhood system. And the overall probability of whole ran¬ 
dom field can be indicated by conditional probability distribution 
describing its local characters. When the noise of MRF 
neighborhood meets the Gauss distribution, a linear equation can 
be obtained which is expressed by a group of airspace element 
gray level, called Gauss-Markov Random Field model. 

Gauss-Markov Random Field (Gauss-MRF) is a steady 
auto-regression course. The covariance matrix of Gauss-MRF is 
positively determined and the neighborhood system is symmet¬ 
rical, with the parameters of the symmetrical neighborhood being 
equal. Expressing texture with Gauss-MRF, namely the gray 
level y(s) of any point s in the picture is the function of 
neighborhood gray level in all direction, which can be shown in 
following condition probability form: 

p(y(s)\all: y(s + r), re N) (2) 

Where, N is a symmetrical neighborhood with s (s is not in¬ 
cluded) being its center and r being its radius. The relation be¬ 
tween the order of Gauss-MRF and the neighborhood system is 
shown in Figure 1. 
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Fig. 1 The model of Gauss-Markov Random Field 

Let S be the point set on M*M net, and 
S = {(i, j), 1 < i, j < M} , the texture y(s), seS, S={(i, j), 1 < i, j 
<M| is a Gauss random process of zero mean value. Then 
Gauss-MRF can be expressed with a linear equation including 
several unknown parameters. 

y(s) = X 9 r [y(s + r) + y(s - r)l + e(s). . (3) 

reN s 


ters Estimation 

If the 2 nd order Gauss-MRF model is adopted, 

I \ T 

Q s =| jw, + y s - n > ■ ■ -,y s+ri + y s -r 4 \ O) 

{rl, r2, r3, r4} = {(0,1), (1,0), (1,1), (1,-1)}, then 0 is a 
four-dimensional vector, and 9 = (0 X , 0 2 , , 9 4 Y ■ 

If the 3rd order Gauss-MRF model is adopted, 

I \ T 

Q s =|^ +ri +y s - v ~-,y s +r 6 +y s -r 6 \ ( § ) 

{rl, r2, r3, r4, r5, r 6 } = {(0,1), (1,0), (1,1), (1,-1), (0,2), (2,0)}, 
then 9 is a ( 6 D) six-dimensional vector, and 

9 = (# 1 , 9 2 , # 3 ,9 4 ,9 5 ,9 6 ) . 

If the 4th order Gauss-MRF model is adopted, 

I \ T 

Qs =|>Wl +T s -, i >• • ->y s +r w + Js-r 10 (9) 


Where, N$is the Gauss-MRFneighborhood system of s, r the 
radius of neighborhood, 0 r the coefficient, and the gauss 
noises array with a zero mean value. Since N is symmetrical, and 
Of = 0_ r , Equation (3) can be expressed as 

y(s) = X e r [ JT (s + '')] + e(s) (4) 

reN s 

Where, (s + r) is a point in a closed annular area S, 
fori' = (/, j), r = (k, l) , we have 

y 1 (s + r) = 

yfs + r) if s + reS (^ 

y[(i + k-l)mod(M + l),(j + l-l)mod(M + l]// s + r«S 


{rl, r2, r3, r4, r5, r6, rl, r 8 , r9, rl0}={(0,l), (1,0), (1,1), (1,-1), 
(0,2), (2,0), (1,2), (-1,2), (2,1), (-2,1)}, then 0 is a (10D) 
ten-dimensional vector, and 

9 = ($[, 92,9s,, $4, 0$, 9§, 9-j, 9%, 9q ,$10) . 

If the 5 th order Gauss-MRFmodel is adopted, 

\ T 

Qs= Ijw, + y s - ri ,■■■, >w I2 + jw 12 \ 0 o) 

{rl, r2, r3, r4, r5, r6, rl, r 8 , r9, rlO, rll, rl2} = {(0,1), (1,0), (1,1), 
( 1 ,- 1 ), ( 0 , 2 ), ( 2 , 0 ), ( 1 , 2 ), (- 1 , 2 ), ( 2 , 1 ), (- 2 , 1 ), (- 2 , 2 ), ( 2 , 2 )}, then 
6 ? is a (12D) twelfth-dimensional vector, and 

9 = ($15 @2 •> @3 > @4 j ^5 j ^6 j ? ^8 j 9g j ^10 j 9j 1 > $12) 

Eqs. (11) and (12) can be solved by Estimating method of 
Least Mean Square as follows: 


If Equation (5) is applied to every point in zone S, the M’ 
equations with respect to {e(s)} and {y(s)} can be obtained 
(Gong 1999; Sheng 2000). 

y( 14) = X^-kt ((1>1) + r ) + e(U) ,., 

reN s 

y(l,M)= Y J 0 l .y l ((l,M) + r) + e(l,M),..., 

reN s 

y(M, 1) = X d r y\ + r) + e{M, 1) 

reN s 

y{M,M) = Y J d r y x m,M) + r) + e{M,M) 

reN s 

Expressed in the form of matrix, it becomes 
y = Q T 9 + e (6) 

Equation (6) is a linear model of Gauss-MRF, of which Q T is 

the matrix about y l (s + f) , and 9 is the feature parameter 
vector of the model to be determined. 

Gauss-Markov Random Field Model and Parame¬ 


e = (Y.QsQl)~'('LQsys) (ii) 

seS^ SGSy 

o=\Y J (y s - 0 T Qs ) 2 ( 12 ) 

N s eS t 

where 9 is the consistent estimation of Gauss-MRF model 
parameters (Gong 1999; Sheng 2000); <7 is the square error of 
parameter estimation. 

When the order of neighborhood is low, there is a certain limi¬ 
tation to describe the feature of complicated pictures with 

Gauss-MRF. 0 can describe more complicated textures with 
high order Gauss-MRF. 

Feature parameter collection and classification 

Feature parameter collection and analysis 

We chose 300 pictures (1:1) of radial texture and tangential 
texture samples of 4 common woods, Pinus koraiensis, Larix 
gmelinii, Fraxinus mandshurica and Quercus mongolica. The 
size of the sample is 12 cm xl2 cm, the resolution ratio is 512 x 
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512 (Fig. 2). Gauss-MRF parameter estimation procedure was 
programmed with Matlab6.5, 2-5 order Gauss-MRF parameters 
were calculated. 



Fig. 2 The pictures of radial texture and tangential texture of four 
wood species 

a is of radial-texture of P .koraiensis; b is tangential-texture of P. koraiensis ; c 
is radial texture of L . gmelinii ; d is of tangential-texture of L.gmelinii; e is 
radial-texture of F. mandshurica ; f is tangential-texture of F. mandshurica ; g 
is radial-texture of Q. mongolica ; h is tangential-texture of Q. mongolica. 


In this paper, as space is limited, only the 2 nd order 
Gauss-MRF model feature parameters B1-B4 of 24 wood sam¬ 
ples are listed in Table 1 and Table 2. Parameter, Bl, B2, B3, and 
B4, describe the texture feature in the direction of 0°, 45°, 90°, 
and 135°, respectively. 

From Table 1 and 2, we can educe the following regulation of 
wood texture. 1) For a certain concrete sample, the maximum 
texture parameter’s gathering direction is main texture direction; 
2) For the samples whose main directions are the same, the finer 
and the clearer the texture is, the greater its corresponding pa¬ 
rameter is, and the smaller the other parameters are. Since the 
sample direction of texture is mainly around 0°, B1 is the great¬ 
est; Bl, B2 of tangential texture are smaller than that of radial 
texture, B3, B4 of tangential texture are greater than that of ra¬ 
dial texture, respectively. Therefore it can be used as the basis to 
classify radial-texture and tangential-texture. In practice, how¬ 
ever, a threshold value is necessary to be determined on the basis 
of a large number of samples. 

In order to determine the validity of Gauss-MRFparameters of 
every order for describing wood surface texture feature, 300 
samples was artificially divided into six texture types, 50 for 
each, and the mean value and square average of the correspond¬ 
ing parameters were calculated, as is shown in Figure 3. It 
showed that the Gauss-MRF parameters of different orders were 
marked different in different texture types. It is feasible to de¬ 
scribe wood texture with these feature parameters. Obviously, the 
higher the order is, the more the parameters are, and the more 
particular the texture is described, but the operation amount will 
be also increased correspondingly. 


Tablet. The second rank Gauss-MRFparameters of radial texture 


Sample 

Bl (0°) 

B2 (90°) 

B3 (45°) 

B4 (135°) 

1 . P .koraiensis 

0.4732584 

0.2470993 

-0.10378 

-0.116577 

2. P. koraiensis 

0.4803995 

0.2515258 

-0.129619 

-0.102304 

3 .P. koraiensis 

0.4788409 

0.2559477 

-0.110502 

-0.124283 

4. L. gmelinii 

0.4995179 

0.3231366 

-0.170195 

-0.152455 

5 .L. gmelinii 

0.4933195 

0.3467532 

-0.165938 

-0.174124 

6.L. gmelinii 

0.4938472 

0.3614616 

-0.203164 

-0.152134 

7. F. mandshurica 

0.4979165 

0.353452 

-0.190948 

-0.160427 

8. F. mandshurica 

0.4955053 

0.3663765 

-0.162626 

-0.199276 

9. F. mandshurica 

0.4994605 

0.3470714 

-0.159778 

-0.186764 

10.(7. mongolica 

0.4989654 

0.3707962 

-0.183878 

-0.185881 

11 . Q. mongolica 

0.499075 

0.3343957 

-0.172211 

-0.161274 

12 . Q. mongolica 

0.497424 

0.3631751 

-0.191612 

-0.168995 


Table 2. The second rank Gauss-MRFparameters of tangential tex¬ 
ture 


Sample 

Bl (0°) 

B2 (90°) 

B3 (45°) 

B4 (135°) 

1 . P. koraiensis 

0.4588136 

0.2509415 

-0.103155 

-0.106618 

2. P koraiensis 

0.4860764 

0.2039888 

-0.1158 

-0.074266 

3. P koraiensis 

0.4738804 

0.1818294 

-0.075555 

-0.080157 

4. L.gmelinii 

0.4928261 

0.2544015 

-0.136343 

-0.110896 

5. L. gmelinii 

0.498788 

0.2733763 

-0.131518 

-0.140647 

6. L. gmelinii 

0.4844843 

0.3396408 

-0.147242 

-0.176888 

7. F. mandshurica 

0.4939181 

0.2921971 

-0.124758 

-0.1614 

8. F. mandshurica 

0.4978635 

0.3419207 

-0.199438 

-0.140378 

9. F. mandshurica 

0.4985287 

0.3194452 

-0.142495 

-0.175541 

10. Q. mongolica 

0.4989141 

0.2748764 

-0.148427 

-0.12537 

11 .Q. mongolica 

0.4951261 

0.2662392 

-0.095775 

-0.165599 

12 .Q. mongolica 

0.4977932 

0.2638029 

-0.156383 

-0.105217 


Feature parameters collection and the texture clustering 

Although we had already calculated Gauss-MRF parameters 
of wood texture of 2 nd to 5th order, which order’s model pa¬ 
rameters are best is still a question. From view point of pat- 
tern-recognition, we must use the separation criterion of texture 
classification to detennine the feature vector of texture rather 
than only according to the difference of a single feature parame¬ 
ter variable (Bian et al. 2000). The classification separating cret- 
erion is defined as 


J d (*) = 


It*. 


V(x! 


,U Is 


(13) 


7=1 


k=1 /=1 


Where C is the number of species, n , and n the number of 
samples in a species, P, and P, the priori probability of the cor¬ 
responding species, and and x\^ the parameters of a 
certain sample in species coi and coj. Let 



k=1 


C 

m = ^ P l m l (15) 

;=l 
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Fig. 3 The 2-5 orders of Gauss-MRFtexture parameters’ average and square difference of six wood texture species 
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C 

S h = Y J P i (m i -m) T {m i -m) (16) 

i =1 

Sw = Y j P i j-f^ (x[°-w,)(4° -my (17) 

1=1 ' k =1 

Where m i is a mean vector of all feature parameter in a 
species, m the population mean vector of all feature parameters 
for all species, S b the parameter distribution matrix of all spe¬ 
cies, and Sw the divergence matrix among all species. Then 

J d =tr(Sw+S h ). (18) 

In Equation (18), the greater the J, is , the easier the species 

is separated. The value of separation criterion is shown in Table 
3 based on Equation (19). 


Table 3. Separation criterion of 2-5 order Gauss-MRFtexture pa- 
rameters of six wood texture species_ 


Orders of 

2nd order 

3rd order 

4th order 

5th order 

Gauss-MRF 





separation criterion 

0.0058 

0.0065 

0.0115 

0.0121 


Table 3 proves that the separation criterion value of 5th order 
Gauss-MRF parameters is the maximum one in six kinds of 
wood texture, so, the 5 th order Gauss-MRFparameters are the 
best one to be used as the feature vector for the texture classifi¬ 
cation. 

In order to verify the validity of the texture parameters in clas¬ 
sification, Hamming competition neural network is adopted to 
classify 300 wood texture pictures with B1-B12 of the 5 th order 
Gauss-MRF parameters, the overall ratio of correctness is 88%, 
which can basically satisfy the requirement for automatic com¬ 
puter classification . 


Conclusion 

The research shows that, by using the Gauss-MRF model, 1) 
the description of the wood surface texture becomes the estima¬ 
tion of texture parameters; 2) the feature of wood surface texture 
described by Gauss-MRF model has a clear physics meaning, the 
gathering direction of the maximum texture parameter is the 
main direction of textures; 3) for those samples with the same 
main direction, the finer and the clearer of the texture is, the 
greater the corresponding parameters are, and the smaller the 
other parameters will be. For the 2 nd order Gauss-MRF, parame¬ 
ter B1,B2 of tangential texture are smaller than that of radial 
texture, while B3 and B4 of tangential texture are greater than 
that of radial texture, respectively. Therefore it can be used as the 
basis to classify radial texture and tangential texture. As for the 
automatic classification of wood textures, the 5 th Gauss-MRFpa- 
rameters were used as the feature vectors in the classification, 
with correctness ratio being 88% as a whole, well satisfying the 
requirement for automatic computer classification. 
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